Theory of the Anderson impurity model: 
The Schrieffer— Wolff transformation re— examined 



Stefan K. Kehrein^ and Andreas Mielke^ 

Institut fiir Theoretische Physik, 
Ruprecht-Karls-Universitat, 
D 69120 Heidelberg, F.R. Germany 



October 5, 1995 



Abstract 

We apply the method of infinitesimal unitary transformations recently introduced by Wegner [1] 
to the Anderson single impurity model. It is demonstrated that this method provides a good ap- 
proximation scheme for all values of the on-site interaction [/, it becomes exact for U = 0. We are 
able to treat an arbitrary density of states, the only restriction being that the hybridization should 
not be the largest parameter in the system. Our approach constitutes a consistent framework to 
derive various results usually obtained by either perturbative renormalization in an expansion in 
the hybridization F, Anderson's "poor man's" scaling approach or the Schrieffer- Wolff unitary 
transformation. In contrast to the Schrieffer- Wolff result we find the correct high-energy cutoff 
and avoid singularities in the induced couplings. An important characteristic of our method as 
compared to the "poor man's" scaling approach is that we continuously decouple modes from 
the impurity that have a large energy difference from the impurity orbital energies. In the usual 
scaling approach this criterion is provided by the energy difference from the Fermi surface. 
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1 Introduction 



In the past few years there has been renewed interest in the Anderson impurity model. This 
model was originally proposed by Anderson [2], for a recent review see [3]. It has been introduced 
to study the well-known Kondo problem, the behaviour of a single magnetic impurity coupled to 
a conduction band of electrons. The Hamiltonian contains the electron band, the energy of the 
impurity orbital together with a repulsive interaction on the impurity site, and a hybridization 
between the band states and the impurity state, 

H = J2 ^kcla^k,a + ^ddUa + J2 ^kicl^a^^a + d^Ck^a) + Ud\dld-d+. (1.1) 
k,a k,a 

In the case of a linear dispersion relation for the band and Vk = V = const, the model was solved 
using a Bethe-ansatz [4]. But if one wants to study the Anderson impurity model in a more 
general situation, one needs a different approach. There are several methods available, most of 
them are reviewed in [3]. The most prominent method amongst them is probably the numerical 
renormalization group developed by Wilson [5] for the original Kondo Hamiltonian and applied 
to the Anderson impurity model by Krishna- murthy et al [6]. 

Recently we applied a new technique to this model [7]. We used continuous unitary trans- 
formations in a form introduced by Wegner [1] to diagonalize the Hamiltonian approximately. 
A continuous unitary transformation yields flow equations for the Hamiltonian, or equivalently 
flow equations for the coupling constants. The approximation we used neglected some additional 
couplings generated by the flow that were not present in the initial model. Unfortunately we 
were not able to obtain quantitative results. The reason is that the continuous transformations 
yields coupled, nonlinear differential equations for the different parameters in the Hamiltonian, 
which we were not able to treat analytically. 

In a more recent work we applied continuous unitary transformation to the well known spin- 
boson model [8]. For this model a unitary transformation exists, the polaron transformation, 
which has been used to treat the Hamiltonian of the spin-boson model approximately. In Ref. [8] 
we used a simple modification of the previous ansatz by Wegner [1] for the generator of the 
continuous unitary transformation. Thereby we were able to construct a continuous polaron 
transformation. This new transformation does not have the disadvantages of the usual polaron 
transformation as it treats the slow bosonic modes in a satisfactory way. In addition, we were able 
to reduce the set of differential equations to a single, non-linear differential equation. This finally 
allowed us to obtain quantitative results which are in good agreement with results obtained by 
other methods. 

In the case of the Anderson impurity model, a unitary transformation similar to the polaron 
transformation is known, the Schrieffer- Wolff transformation [9]. It has been introduced to 
eliminate the hybridization between the electronic states in the band and in the impurity orbital 
present in the Anderson model. Thereby it renormalizes the impurity energy and the repulsive 
interaction. Furthermore it generates a spin-spin interaction between the impurity electron 
and the conduction band electrons. This interaction is responsible for the Kondo effect. In 
this treatment one usually takes only the terms into account which are of second order in the 
hybridization. Other interactions are generated as well, but they are neglected. 

The Schrieffer- Wolff transformation has some disadvantages. First it is equivalent to a second 
order perturbational treatment of the hybridization term in the Hamiltonian. Therefore the 
general validity of the result is unclear. In particular the Anderson impurity model in the Kondo 
regime is mapped onto a Kondo problem with an effective band width of order of the conduction 
band width. This is known to be wrong since the high-energy cutoff in the Anderson impurity 
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model cannot be larger than of order the on-site interaction U [10]. The second problem with 
the Schrieffer- Wolff transformation is that energy denominators occur, which become zero if the 
energy of the impurity level lies in the conduction band. If this is the case, the Schrieffer- Wolff 
transformation is ill-defined. This problem is similar to the problem of the treatment of the 
slow bosonic modes in the polaron transformation. Therefore we expect that a modification of 
our old treatment of the Anderson impurity model should be useful. In the present paper we 
construct a continuous, or infinitesimal Schrieffer- Wolff transformation in order to eliminate the 
hybridization terms in the Hamiltonian. We will show how such a continuous Schrieffer-Wolff 
transformation can be constructed systematically. Our old approach can be modified so that the 
number of additional interactions generated by the unitary transformation is reduced. 

It is clear that applying infinitesimal unitary transformations to a given Hamiltonian is a 
non-perturbative method. The advantage of our present approach is that our method treats 
the Anderson impurity model consistently within one framework independent of whether the 
on-site interaction U or the hybridization F is the larger parameter. We are able to reduce the 
problem to two coupled non-linear differential equations for the impurity orbital energy and 
the on-site interaction. These are solved approximately yielding self-consistency conditions for 
these two quantities. Finally the antiferromagetic spin-spin interaction can be calculated. It 
is demonstrated how standard results from e.g. renormalization theory can be obtained in this 
conceptually new framework. Whereas other methods are not applicable in the whole parameter 
space or need additional assumptions, continuous unitary transformations are conceptually sim- 
ple and no physically relevant restrictions or additional assumptions are needed. One can hope 
that the flow equations approach will be useful too for other problems with less well-established 
results. In so far it is important to study the method using a well-known problem. 

Independent of Wegner [1], Glazek and Wilson [11, 12] have recently also proposed to use 
continuous unitary transformations to construct renormalization group equations for effective 
Hamiltonians in quantum field theory. Wilson et al. [13] applied this method to quantum 
chromodynamics. Although the general idea is similar, there are some differences between the 
approach of Glazek and Wilson and ours. Their goal is to eliminate "far-off-diagonal" matrix 
elements in a given Hamiltonian, which means off-diagonal matrix elements connecting states 
that are energetically far from each other. This means that the final Hamiltonian has a banded 
structure. In contrast our goal is to eliminate matrix elements such that the final Hamiltonian 
is diagonalized approximately or block diagonal. 

Although it is possible to construct flow equations such that the Hamiltonian becomes diag- 
onal, it is not possible to solve these equations and to calculate the eigenenergies. Therefore we 
use the continuous unitary transformation in order to eliminate some of the matrix elements. In 
our case the hybridization is eliminated and the flnal Hamiltonian is block diagonal. The final 
Hamiltonian does not contain matrix elements connecting states with a singly occupied impurity 
orbital to states for which the impurity orbital is either not occupied or doubly occupied. If now 
the impurity site is occupied with no or two electrons, the spin on the impurity site is zero and 
the additional antiferromagnetic interaction vanishes. In these cases the problem is essentially 
solved with respect to static properties. In the regime where the impurity site is singly occupied, 
the problem is reduced to a usual Kondo problem, which may be solved by various methods. 

Our paper is organized as follows. In the following section we derive the general flow equations 
for the coupling constants of the Anderson impurity model. In section 3 we illustrate the method 
in the case of a vanishing interaction of the electrons on the impurity. It is shown that in this case 
the flow equations yield the exact solution. Furthermore we introduce an approximation which 
still gives the exact result for the case of vanishing interaction and which is applied in section 4 
to the case of a non-vanishing interaction. We calculate implicit equations for the renormalized 
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interaction and the renormalized impurity energy. These equations can in principle be solved 
for any given density of states and hybridization. Our results are compared with the results 
obtained by Schrieffer and Wolff [9]. In section 5 the antiferromagnetic spin-spin interaction is 
calculated. It can be used to determine e.g. the Kondo temperature. In section 6 we discuss the 
method of continuous unitary transformations as compared to a single unitary transformation 
in the Schrieffer- Wolff paper. Section 7 contains a discussion of our method as compared to the 
well-known "poor man's" scaling approach. The last section contains the conclusions and an 
outlook on other problems. 

2 The flow equations for the Anderson impurity model 

Our starting point is the Hamiltonian for the single impurity Anderson model (1.1), which we 
write in a normal ordered form 

H = J2^k- A,a('k,a ■■ + Yl ^ddld„ + J2 MA,ada + d^Ck^a) + Ud\dld^d+ + Eq. (2.1) 
k,a <^ k,a 

We introduced a normal ordering for the band electron operators ^ and c;t,(T- It is defined by 
• '^'ka^k.ij '■= c\ ^Ck,cr — where Uk = (exp(/3efc) + 1)^^ is the occupation number of the band 
state with wave vector k. We let the Fermi energy equal to zero. The reason for the normal 
ordering is that additional interactions, which will be generated by our procedure and which we 
neglect, should be written down in a normal ordered form as well since otherwise the ground 
state expectation value of such additional contributions does not vanish. A detailed discussion 
has been given in [7]. In contrast to [7] we do not introduce a normal ordering for the impurity 
electron operators dj. and d^r- The reason is that in our approximation contributions containing 
such operators are not neglected if they have a non-vanishing expectation value in the ground 
state. It is possible to introduce a normal ordering on the impurity site as well, but the results 
are not changed. Due to the normal ordering a constant Eq = 2 e^rifc has been introduced in 
the Hamiltonian. 

We now want to apply a general continuous unitary transformation to the Hamiltonian. Such 
a transformation is defined by a generator r] that depends on a continuous variable, which we 
call l. The continuous unitary transformation is defined by ^ = [7]{i) , H with the initial 
condition i?(0) = H in (2.1). In order to simplify notation we will denote the initial values of 
and U by 

'i? e^{l = 0), = U{£ = 0) (2.2) 

and the asymptotic ( "renormalized" ) values ioi £ ^ oo by 

er=l^,(^ = oo), U^'^'ui£ = oc). (2.3) 

If these parameters appear without an argument this will imply that they are to be considered 
as functions of £. 

We choose r) to be of the form 

V = Y.^k{ci^ada - d\ck,a) + ^ %,g(4,aCg,^ - c\^„Ck,a) 
k,a k,q,a 

+ E ^k^ {clJlad-ada - dldl^d^^Ck,.) (2.4) 
k,a 
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The commutator of rj and H is easily calculated, we obtain 

k,q,a 

k,a k,a 

+ Uj2vk {dtdl^d^aCk,a + ci^a^-ad-ada) 

k,a 

- J2 ^k,g{ek - eg){: c^^Cg^a : + : claCk,a 

k,q,a 

+ 2 ^ r]k,gVg{dlck,a + 4,a'^'^) 
k,g,a 

- J2 Vk^ i^k - £d) {dldl„d^„ck,a + c\jl^d^M 

k,a 

- 2j2v''k^Vkdidl,d^,d, 

k,a 

+ J2 ^k^^li- (^k,ad-ad-aCg,a ■ + : C+ ^^C?L^C?_^Cfe,^ : - : cl^^dl^d^Cg^^^ : 
k,g,a 

~ '■ '^g,a'^-a'^cr^k,-a '■ ~^k^a^g,-a'^-cr'^cr ~ dljd^_^Cg^-crCk,a) 

+ 2j2v'i^Vknkdl,d^^ 

k,a 

+ Uj2r)^k\dU'-ad-aCk,a + clJlad.M- (2.5) 

k,a 

Many additional couplings are generated which did not occur in the original Hamiltonian in 
(2.1). But some of these terms can be eliminated by a suitable choice of r]. Let us first consider 
terms containing : ^Cg^a^ + c'^^^^k^a '■■ Such terms do not occur if we choose 

Vk,q{^k - eg) = ^iVkVg + rigVk). (2.6) 

Similarly, terms containing operators of the type d^^dL^d^fjC^^^ + c^. ^d^_^d^fjdfj do not occur if 
we choose 

Uvk + {ea - ek)v?^ + Uri^ = (2.7) 

(2) 

These equations may be used to determine and . The only contribution that does occur 
additionally is the term 

Yl ^k^^li' A,ad-ad-aCg,a : + : c\^„dl„d-„Ck,a ■ " : cl„dl„d„Cg,-„ : 
k,q,a 

- ■■ 4,<Td-adaCk,-a ■ -cl^^c\_„d-„d„ - d^dl^Cg-^Ck^a) ■ (2-8) 

Notice that the contributions containing : c^, ^Cg^fj + c''g^(j^k,a '■ or • d]jd^_^d^fjCk^i^ + ^d^_^d^fjdfj : 
have been classified irrelevant in our former approach [7], whereas the term in (2.8) is marginal 
in some of the fixed points [7]. In principle it has to be included in the Hamilonian in (2.1). 
But the commutator of this term with r] does not yield contributions to the other terms in the 
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Hamiltonian. Therefore we do not take it into account in our first analysis of the problem. 
But it is clear that this additional term is important. A part of it yields the antiferromagnetic 
interaction between the impurity spin and the spins of the band electrons that is responsible for 
the Kondo effect. We will come back to this interaction later. Our first goal is to calculate the 
flow equations for the parameters in the Hamiltonian (2.1). Prom (2.5) we obtain 

^ =r]k{ed-(^k) + '^Y.^k,pVp (2.9) 

^ = 2%F, (2.10) 
^ = -2 5]%F, + 2^4'V,n, (2.11) 

= -4 5]4^V, (2.12) 



dl 



k 

dEp 
dl 



4 5]%Ffcnfc (2.13) 



The last equation yields directly Eq = 2 Y,k ^-k'^k^ which is the energy of the filled Fermi sea. In 
the following we are interested in the thermodynamic limit. For large iV, the number of states in 
the band, one has Vk oc N'^l"^. Thus, r]k must as well be of the order N'^l"^, and the derivative 
of e/j with respect to I is of order N^^. For large values of N the band energies do not depend 
on I. This should have been expected. The thermodynamic bath of electrons is not affected by 
the single impurity. This means that the global density of states in the band is fixed. But this 
does not mean that the local density of states is fixed as well. In contrary, one should expect 
that the local density of states near the impurity site is affected by the impurity. We will come 
back to this point in the discussion. 



3 Vanishing interaction [7 — 

To illustrate the advantages of our method, let us first study the case C/ = 0. Then we have a 
quadratic Hamiltonian that can be solved exactly, see for example Ref. [14]. We will show that 
our method yields the exact solution in this case. For [/ = the flow equations simplify to 

^ =%(ed-eA;) + 2^%,pV^p (3.1) 

^ = -2E%^^ (3.2) 

The flow equations are exact in the case [/ = since the neglected terms in (2.5) vanish in this 
limit. We let 

r]k = Vkf{ekA) (3.3) 
and introduce 

J{e,l) = Y,V^5{e-ek). (3.4) 
k 
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In the literature one often introduces the parameter 



T = 7rp{eF)VkA0f = 7rJ{eF,0), (3.5) 

where p{eF) is the density of states at the Fermi surface. The flow equations for and J{e,l) 
are 

ded 
dl 



-2 j def{e,i)J{e,e) 

J e — e' 



(3.6) 



d£ 



(3.7) 



For the integral in the last equation one has to take its principal value. The set of flow equations 
may be solved if one introduces a function 



G{eJ) = j de' 



J{e',i) 



e-e' + G{e, 



(3.8) 



Taking the derivative with respect to i, we obtain an implicit equation for this derivative, which 
can be solved. The final result is 



dG{e,. 



1 + 



J(e',. 



+ Gi€JW 



-1 



de f{e , i)J{e' , 



eg-e' + G(e, 
e-e' + G{e,. 



Calculating the derivative of G{e,i) with respect to e we obtain similarly 



1 + 



dG{e,i) 
de 

This yields 

dG{e,i) 
di 

1 + 



1 + 



J{<'J) 



dG{e,l) 
de 



2 / de'f{e\i)J{e' 



{e-e' + G{e,i))\ 

ea-e' + G{eJ) 



e-e' + G{e,l) 



Comparing the right hand side with the derivative of with respect to ^, we obtain 



ded 
d£ 



dG{e,£) 

at 



1 + 



dG{e,t) 
de 



This equation can be integrated and the final result is 
ed{i)+Giedii),l)=el 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



In the last step we used G(e, oo) = 0, which follows directly from J(e, oo) = and holds for an 
appropriate choice of f{e,£). Solving for erf(^), we obtain 



ed{i) 



I 



de 



(3.14) 



Again, for the integral on the right hand side we have to take its principal value. This means 
that we have to choose J{e,l) and therefore f{e,i) such that the principal value exists for all £. 
We obtain the value of e^, if we let ^ = and solve for e^. 



J(e,0) 



As a simple example we take a semi-circle 
2F2 



VD^ - e^ 



(3.15) 
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where 2D is the band width and V = ^J2k ■ Here we have F = 2V'^ /D. The main reason for 
this choice of the hybridization is that all the integrals can be worked out in closed form in the 
sequel. But it should be noted that it is a main advantage of our approach that it can be used 
for arbitrary functions J(e, 0), in particular for any distribution of the density of states in the 
conduction band. However, if one chooses a linear dispersion relation and constant hybridization 
Vjfc, that is J(e, 0) = |^©(D — |e|) as usually done in renormalization group treatment of the 
Anderson impurity model, one must be careful due to the discontinuous behaviour close to the 
band edge. The self-consistency equations in the flow equations approach will then generally 
have more than one solution, however, the actual solution of the differential equations chooses 
the correct one. Close to the band edge one expects unphysical behaviour anyway due to the 
unpliysical choice of J(e, 0) and neither approach should be trusted. Therefore it is natural in our 
approach to choose a function J(e, 0) that is continuous at the band edge as should be expected 
on physical grounds anyway. 

Let us come back to our example introduced in Eq. (3.15). We have to distinguish between 
the case where lies in the band and the case where it lies outside the band. We first consider 
the latter case. The integral can be calculated and we obtain 

^ 6^ + § (sign(ef -D^- ) = (3.16) 



d 



if 



-R 



> D. This equation yields a simple quadratic equation for e^, which has always two 
solutions. If r < D at most one of these solutions lies outside the band, li T > D and 
> V — D, there is a single solution for outside the band, but if < F — D we obtain 
two solutions outside the conduction band. 

The situation is much simpler when lies inside the band. The integral in (3.14) has to be 
interpreted as its principal value and we obtain 

4-^d-^4 = ^- (3-17) 
The only solution is 

4 = (3.18) 



1-75 



lies inside the band if F < D — 
The fact that for a sufficiently 



The various cases are shown in Fig. 1. 

arge value of F (F > D in our example) two solutions for 
exist, is generic. It holds for any J(e, 0) with a connected support of length 2D. It is clear that 
with the present approach of flow equations only a single solution can be obtained. Nevertheless, 
the second solution is of physical importance. It is possible that a localized state develops from 
the original band states that has an energy which lies outside the band. Such a state cannot 
be obtained within the present formulation of the flow equations. In the case f/ = one can 
introduce a different representation of the Hamiltonian and of r] that includes a localized band 
state explicitly. Since we are at present not able to deal with similar problems in the case ?7 > 0, 
we restrict ourselves to the parameter regime F < D (in fact later we will need F < D 12). This 
is reasonable from a physical point of view since we do not expect that the hybridization is the 
largest parameter in the system. 

Eq. (3.14) is obtained as well if one uses two simple approximations to the flow equation. 
The first approximation neglects the terms proportional to c\ ^Cq^a in the Hamiltonian that are 
generated by the transformation and consequently one neglects such terms in -q as well. This 
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is an approximation that can be justified from a physical point of view, since these terms are 
irrelevant in all fixed points [7]. Then the second term in (2.9) vanishes and the equation for 
is linear in Vk- Similarly, the second term in (3.7) vanishes, whereas (3.6) remains unchanged. 
Both equations together yield 

^ = - / de^. (3.19) 
a£ J Ed — e 

Furthermore we assume that erf(^) converges rapidly to ef, so that we can replace ed,{£) with 
on the right hand side (3.19). This yields (3.14). We will use similar approximations for [/ > 
as well. Although it is possible to choose r/ in such a way that only very few new terms are 
generated, the flow equations become very complicated. In order to be able to analyse the flow 
equations, one has to neglect higher interactions. This is often possible due to physical reasons. 
One would like to understand why the self-consistency condition obtained by replacing by 
on the right hand side of (3.19) yields a good approximation to the exact solution of (3.19). 
To discuss this point let us introduce a special choice of /(e,£). Since we want J(e,£) to vanish 
in the limit ^ — > oo, a natural choice would be f{e,i) = — (e^ — e). But for finite U we will 
have to make a different choice for /(e,^) in the next section. For consistency we therefore take 
f{e,tj = —(erf — e)'^/(4e^). This obviously works as well and it is easy to see that in the present 
case both choices are essentially equivalent. 

The following argument holds in both cases. Unless e = e^, J(e, 1) decays exponentially on 
a scale set by ^ oc ^\l{^d — e)^- If lies outside the band it will tend to exponentially and 
the approximation « on the right hand side of (3.19) is justified. On the other hand, if 
lies inside the band, we can estimate the relevant ^-scale on which changes by calculating the 
ratio of the total change of to its derivative with respect to I for small i. This shows that 
changes on a scale set by ^ oc e^/D^, i.e. much faster than J{e,£) for values of e near the Fermi 
energy. Therefore can be replaced by its renormalized value on the right hand side of (3.19). 
We will use the same approximation in the next section to discuss the case ?7 > 0, it can be 
justified in the same manner. 



4 Non— vanishing interaction U > 

With the approximations introduced at the end of the last section, the flow equations for [/ > 

may be written in the form 

^=%(ed-efc), (4.1) 
^ = -2 5:%F, + 2^4^V,n,, (4.2) 

k k 

^^-4i:.fV.. (4.3) 

k 

According to (2.7) we take r]^^'^ = —Ur]k{ed ~^k + U)~^ and as above r/^ = Vkf{ek,£)- We assume 
that J2k th^aA, the renormalized is unique for [/ = 0. We expect that it is unique 

for ?7 > as well. With these assumptions we proceed as in the previous section. We introduce 
J{e,£) as in (3.4) and obtain the flow equations 

^:l^=2f{e,i)J{e,i){ed-e), (4.4) 
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ded ^ r dJ{e,i)ed-e+{l+n{e))U 

d£ J ^ d£ {ed-e){ed-e + U) ' ^ 

dU ^ f dJ{e,£) U 

d£ J d£ {ed-e){ed-e + Uy ^ ' 

n(e) is the Fermi distribution. The first equation may be used to parametrize J{e,£). A suitable 
parametrization is 

J(e, £) = J{e, 0) exp (- f dl\ . (4.7) 

We will see that with this choice the hybridization flows to zero for all e, in particular also for 
e = or e = + U^. The reason is that — decays like as we will see below. J(e^, £) 
decays algebraically to zero. Furthermore, Eq. (4.7) corresponds to the following function f{e,£) 

/(M) = (4.8) 

The reason for this choice of /(e, £) or J(e, £) is that now no pole terms appear in the integrals on 
the right hand sides of (4.5) and (4.6). The denominator in (4.8) is just introduced for convenience 
so that limits like liuiu^;^ can be performed in all the equations without difficulties. Later we 
will come back to the question of other parametrizations /{€,£). In fact Eq. (4.8) belongs to a 
class of parametrizations that all give the same physical results, whereas other parametrizations 
lead to divergencies or J(e, £) does not flow to zero everywhere. 

Let us consider for a moment the simplified case lim[/_^oo- The equation for takes the form 

^ = I deJ{e, 0)(1 + n{e)){ed - e) exp j\ed - efd£'^ . (4.9) 

This equation is very similar to the fiow equation for the renormalized tunneling frequency in 
the spin-boson problem [8]. The asymptotic behaviour can be obtained as in [8], one finds 

~ ^d{() oc for large I if lies inside the band. Otherwise it decays exponentially. (4.9) 
shows that I/a/Z plays the role of an effective band width if becomes smaller than the 

original band width of J(e, 0). The effective band width is reduced with increasing £. This is 
similar to a renormalization group procedure, where modes with high energies are integrated out. 
In our case, these modes are decoupled from the system. The analogies with renormalization 
theory will be worked out in more detail in section 7. 

For finite U the situation is somewhat more complicated, but the results are similar. The 
effective band width is l/y/l and leads to an asymptotic behaviour — U{£) = Cii^'^ and 
ef — f-d{£) = —6*2^^2 for large £ with some constants Ci and C2. This again holds if and 
ed + U lie inside the band. Ci and C2 are positive if lies below and + [/ lies above the 
Fermi energy. One possibility to obtain these results is to make the ansatz U{1) = + Ci£°' 
and ed{£) = ^d ~ ^'2^^- Inserting these expressions in the flow equations one shows easily that 
a = P = ^ is the only possible solution. We now replace U and by their asymptotic values on 
the right hand side of (4.5) and (4.6). Both equations can be integrated and we obtain 

e?-e+(l+n(e))[/" 



ed{l) = e^- [deJ{e,£ f \ ' ""^^1^ ^^-^O) 

\e^-e)ief-e + URy 



m = U- + 2jdeJie^£) (4.11) 
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These equations are good approximations to the solution of the flow equations (4.5) and (4.6). 
They give the correct asymptotic behaviour and a numerical integration of the flow equations 
shows that the true solution differs only slightly from the approximate value. It should be 
noted that this is a different level of approximation than the previous restriction to a certain 
set of interactions included in the flow equations. This restriction was a physical approximation 
whereas the approximate solutions (4.10) and (4.11) can be controlled by solving the original 
differential equations numerically. We have done that too and always found very good agreement. 

For ?7 = we showed that (4.10) yields the exact result. For [/ = oo, (4.10) is correct up 
to terms quadratic in J{e,l). This can be seen if one notices that (4.10) is the exact solution of 
a set of equations similar to (3.6) and (3.7) but with J{e,i) replaced by J(e,^)(l + n(e)). An 
additional argument to justify this approach is similar to the one given at the end of the previous 
section. The relevant ^-scale for changes of and U is smaller than the scale on which J(e, i) 
varies. The crossover to the asymptotic behaviour occurs for £ > D~^, whereas J{e,£) does 
not change too much on this scale. Taking £ = 0. (4.10) and (4.11) yield the self-consistency 
conditions for and U^. Let us mention that the results for ef and obtained from (4.10) 
and (4.11) do not depend on the special choice of f{f-,£) in (4.8). Nevertheless (4.10) and (4.11) 
are only good approximations to the flow equations (4.5) and (4.6) for special choices of f{e,£) 
like the one in (4.8). The important point is that J{e,£) has to be chosen so that the principal 
value of the integrals in (4.10) and (4.11) is well-defined for all values of £. This is clearly 
true for J{e,£) given in (4.7). Our results here do not depend on the details of the continuous 
unitary transformation. But the transformation has to be chosen such that the flow for all the 
parameters in the Hamiltonian is well-defined. In section 7 we will see an example of a different 
parametrization of J(e, £) where this is not the case. 

Let us again consider the case J(e, 0) = ^p^/D^'^^^. The equation for does not contain 
a factor n(e) and the integral is easily evaluated. The result is 

2r 



D 



- eilef + - D)sign{e^ + C/^)^(eJ + C/^)2 - L»2 



(4.12) 



This equation shows that if lies below the Fermi energy and + lies above the Fermi 
energy, then is larger than the initial value . This is also true if both, and + lie 
in the band. It this case we simply obtain 

= (4.13) 

^ D 

On the other hand, if and ef + lie above the energy band, is smaller than . 



Symmetric Anderson model 



In the symmetric case = — 2e^, the fiow equations yield U{£) 
conditions give = —2ef as it should be. Therefore we have 

2r 



D 



-2ed{£) and the above 



(4.14) 



This equation is very similar to the one obtained for f/ = 0. If 



< D — 2T the renormalized 



value is ef = e^J{l - 2r/D) and lies inside the band. If > D - 2r we obtain a quadratic 
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equation for e^. For 2T < D this equation has a single solution outside the band, whereas for 
2r > we can obtain two solutions of the self-consistency equations outside the conduction 
band similar to the case U = 0. 



Asymmetric Anderson model 

In the general case {U / —26^), we have to calculate the integral in (4.10). It contains a factor 
n(e) due to the normal ordering we introduced in the Hamiltonian. Therefore the renormalized 
impurity energy ef depends on the temperature and the chemical potential. We let T = and 
ep = 0, so that n(e) = 1 — ^(e). In this case the integrals can easily be evaluated explicitly. We 
have to distinguish various cases of whether the impurity orbital energies lie inside the conduction 
band or outside. 



R 



> D: 



We obtain 



2D 



2eJ -U^ + sign(e^ + U^^)^ {e^j + U^f - D"^ [1 - - arcsin 



D 



-sign(e^)y^ {(.^Y ~ arcsin 



D 
7r 



(4.15) 



In the limit = we have = and the condition for is the same as in the previous 
section. In the limit U = oo, we obtain a single equation for ej, 



^rl — ^rl + 



r 

2D 



sign{e^)^{e^r-D^ ( 3 - - arcsin ^] + 'lD-3e^ 



2 



(4.16) 



In both expressions, 



.R 



> 



. The impurity orbital energy is pushed away from the band 



as it should have been expected. 



.R 



> D: 



We obtain 



— ^d 



2D 



2e'i-U^ + ^^D^-{4)Hn 



'^D^-{efY+D' 



^R\ 



+sign(ej + U^)^J{e^ + - D"^ ^"^^^^^ 



D 



.R 



(4.17) 



If we let [/ = oo in this case, this expression simplifies to 

^D^-{e^y+D' 



^d — ^d 



V 

2D 



ie^-lD + lP^-ie^fln 



^R\ 



(4.18) 
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Finally for 



4 + - In 



2D 



e e 



R 



(4.19) 



This is shown in Fig. 2. Eq. (4.19) contains a logarithmic singularity on the right hand 
side for 0. If is negative and sufficiently far away from the Fermi energy, is 

negative as well and increases with increasing e^. At some (still negative) value of the 
renormalized impurity energy jumps discontinuously from a given value below the Fermi 
energy to a value above the Fermi energy and then increases further with increasing e^. 
This behaviour can be deduced from the numerical solution of the differential equations 
and is depicted by the full line in Fig. 2. However, it is difficult to obtain reliable numerical 
results in this regime, never reaches the Fermi level except for the trivial case where 
F = and = 0. 



.R 



+ U 



R 



< D: 



We obtain 



F 

2D 



2ef-U^ + -^D^-ief)Hn 



i7ri 



'^D^-ie^ + W^r+D' 



.R 



(4.20) 



This expression contains two logarithmic singularities, it diverges if either ^ or 
ef + 0. As a consequence, both ef and ef + cannot approach the Fermi energy 

as long as F / 0. For e? <^ <^ D one finds in particular 



^d 



4 + 



In 



TT 



u 



R 



,R 



jjR 

+ ro(-). 



(4.21) 



Again the solution of Eq. (4.21) is non-unique for some range of the initial parameter e^. 
In fact Eq. (4.21) is well-known from renormalization theory [6, 10, 16]: In the valence 
fluctuation regime of the asymmetric Anderson model one has to replace by an effective 
impurity orbital energy (we use the notation from [6]). E"^ can be obtained as the 
solution of the following equation 



-Ea{T^) -E, 



(4.22) 



with 



-EdiT) = - 




(4.23) 



One easily checks E'^ = e j with e j from the flow equations approach. This result for e j 
will play an important role for the calculation of the Kondo temperature in this regime in 
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Sect. 5. It will then become apparent why it is important to find the effective impurity 
orbital energy from renormalization theory to be identical to our asymptotic value e^. 

Finally it should be emphasized that the flow equations immediately give the correct high- 
energy cutoff in Eqs. (4.19) and (4.21). The smaller of the two parameters U and D appears 
in the logarithm and this comes about here very naturally. 

These results, especially the self-consistency conditions in the general form 
J , f,.u. -e+(l+n(e))C/^ 



/d..(.o) f,-;'i""':»';; (4.24, 



can be compared with the result of the Schrieffer- Wolff transformation. If one applies a usual 
Schrieffer- Wolff transformation to the Hamiltonian, the impurity energy and the interaction 
U are renormalized as well. But the result is a simple result of a second order perturbational 
treatment. It may be obtained from our self-consistency conditions if they are solved recursively 
to first order in J(e, 0). This simply means that on the right hand sides of (4.24) and (4.25) 
the renormalized values are approximated by the initial values. As already mentioned our result 
is exact if [/ = 0. In this case the self-consistency condition for can be solved iteratively. 
This corresponds to summing up the whole perturbational series. Similarly, the self-consistency 
conditions (4.24) and (4.25) can be solved recursively though this will in general not give the 
exact result. But it is seems that a large part of the perturbational series is summed up when 
one solves these equations due to the same reasons that were already mentioned at the end of 
section 3. 

One of the main results that we have obtained for J(e, 0) = ^^^VD"^ — <? was that the 
renormalized values and + U^^ behave discontinuously at the Fermi energy as a function 
of the initial values e\ and + ?7^. This is a consequence of the fact that due to the normal 
ordering a factor n(e) appears in (4.24). This result is generic and holds for a general fimction 
J(e, 0). A similar effect occurs at the band edge if J(e, 0) is not a continuous function at the 
band edge. Then we obtain singularities in the integral in (4.24) if or + U^^ approach the 
band edge. Consequently these quantities behave discontinuously at the band edge as a function 
of the initial values. 

There is another interesting point that we would like to mention. In the symmetric case 
U = —2ed the renormalized value of U does not depend on the temperature. If the system 
deviates only a bit from the symmetric case, (4.24) and (4.25) show that the system is pushed in 
the direction of the symmetric case. This can be seen if one takes = —U/2 + 5 and expands the 
right hand side of (4.24). The renormalized value of 5 is smaller than the initial value of 5. In this 
sense the symmetric situation is stable. Near the symmetric point the temperature dependence 
of the renormalized values will be weak. But in the general case the renormalized values of 
the impurity energy and the interaction will depend on the temperature. For small enough 
temperature the integral in (4.24) can be evaluated using the usual Sommerfeld expansion. This 
yields 



4 = ^^ 



This shows that for 



> 



.R 

f rl 



we obtain e?(T) > e?(T = 0), whereas for 



< 



we obtain e^(T') < e^{T = 0). Generally Eq. (4.26) is a good approximation only for T < 



.R 



.R 
^d 
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> 



When T becomes larger one can show that e^(T) will decrease as a function of T for 

. The temperature dependence of ef leads to a weak temperature dependence of U^. In 
the special case J(e,0) oc VD^ - e2 we obtain a temperature dependence of only if ef or 

+ lie outside the band. If , + U' 
the temperature. 



< D, (4.12) shows that U does not depend on 



5 The induced spin-spin interaction 

So far the contribution (2.8), which is generated by the unitary transformation, and which 
therefore has to be included in the Hamiltonian, was not taken into account. This interaction 
gives rise to a spin-spin coupling term 

-2E<?(^1^^^«) ■ (5-1) 

k,q 

with 

the so called potential scattering term 

^E<?(4''^«)(^i^^)' (5-3) 

k,q 

and a term 

k,q 

The final Hamiltonian contains no couplings between states that have a singly occupied impurity 
orbital and states for which the impurity orbital is either empty or doubly occupied. Whereas the 
spin-spin coupling (5.1) acts only on the part of the Hilbert space of states with a singly occupied 
impurity orbital, the term (5.4) vanishes on this part of the Hilbert space. It is important if 
the impurity orbital is either empty or doubly occupied. In this sense, these two terms are 
conjugate to each other. In fact one has a simple interpretation for these couplings in the 
symmetric Anderson model. Whereas the asymmetric Anderson model has only the usual SU(2)- 
spin symmetry, the symmetric Anderson model has an additional SU(2)-pseudo-spin symmetry. 
Introducing the wave vector tt, which has all components equal to tt, the symmetric energy band 
has the symmetry ei~ = — ejr-A;, Vk = V^-k- For = —U/2 the Hamiltonian also commutes with 
the operators 



Sz = 1 - d\d+ - did- + E(l - 4,+Cfc,+ - 



k 



S+ = d+d^ + ^c„_k,+Ck- (5.5) 



k 



k. 



These operators form the second SU(2) symmetry mentioned above. The potential scattering 
term (5.3) and the term (5.4) together can be written as a pseudo-spin interaction. It is clear 
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that if the original Hamiltonian has these symmetries, the transformed Hamiltonian has these 
symmetries too. Therefore the term (5.4) is present if the corresponding spin-spin interaction 
is present. Although in the asymmetric case the Hamiltonian does not have the additional 
symmetry, the two terms (5.3) and (5.4) have the same interpretation. The only difference is 

(21 

that now the coupling constant ^ is not symmetric with respect to a transformation k ^ n — k. 

In the remaining part of this section we want to discuss the regime where the impurity 
orbital is singly occupied. Since ipli^ipd = 1 in this regime, the pseudo-spin interaction reduces 
to a scattering of band electrons. Such a term has already been neglected and therefore this 
contribution is not taken into account. 

(2) 

In spite of the fact that the additional couplings ^ are small compared to the other para- 
meters of the Anderson model, they cannot be ignored. In the regime where the impurity orbital 
is singly occupied, ef < ep and e^ + U^ > ^F, even a small antiferromagnetic spin-spin coupling 

(2) 

at the Fermi surface Vj^^^,^ < gives rise to the Kondo effect for low temperatures. In the Kondo 
model the Kondo temperature can be defined as Tk = {2'KXimp{T = 0))^^ [4] where Ximp{T = 0) 
is the impurity contribution to the susceptibility at zero temperature. Let us remark that other 
definitions of the Kondo temperature can be found in the literature, the definition by Wilson [5] 
is somewhat different. For a detailed discussion see e.g. [4]. Based on a Bethe-ansatz solution, 
Tsvelick and Wiegmann argue that the Kondo temperature Tk is given by 

knTK = exp [-$(2p(e^)F,(^;, (5.6) 

with the universal function [5] 

^(y) = A-^ln|?/| + 0(?/). (5.7) 

\y\ ^ 

For the Kondo problem D is the conduction band width. In the Anderson impurity model one 
knows that D in (5.6) has to be replaced by an effective band width Deflf that cannot be larger 
than [10]. If one follows the perturbative calculation of e.g. the susceptibility in the Kondo 
problem, one notices that the breakdown of the perturbation expansion is due to the matrix 

(2) 

elements W '■. These describe the scattering of an electron from the Fermi surface with the 
impurity to some wave vector q and then back to the Fermi surface. For this reason one is not 

(2) (2) 

only interested in the coupling at the Fermi surface V^^\.^ , but also in Vj^^^^ since this determines 
the effective band width of the associated Kondo problem. 

(2) 

Let us now calculate the matrix elements ^ in the flow equations approach. We already 
mentioned that the additional couplings (2.8) do not lead to a contribution in the equations for 

and U. Therefore we calculate the coupling constant simply by integrating the coefficient in 
front of the interaction term in (5.1). To be consistent with the notation in our previous paper 
we symmetrize this coefficient. We have to calculate 

/•OO 

VS= dlU^W, + r^f)V,). (5.8) 





Using (2.7) and (3.3) to replace rfj^ and furthermore (4.4), we obtain 

1 .oo / d\nJ(ek/) d\nj(e„l) , 
-- / diVkV.U — - + — - . (5.9) 



ed - ek){€d - ek + U) (ed - eq){ed - eg + U) 
The ^-dependence of Vjt is obtained from (4.1). Using again (3.3) and (4.4) we obtain 
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This yields 

n^,? = -^^fe(0)F,(0)|^"d£C/(J(e,,^)J(e,,0)J(e„£)J(e„0))-^ 

{^d - ^k){^d - + U) {ed- eq){ed- eg + U) j ' 
Using the parametrization for J{e,i) introduced in (4.7), this expression simphfies to 
y(2) _ lvu(0)V(0) r dni ^^'' ~ ''^^'^ -^k + U) + {ed - eg) (erf -e^ + U) 




Let us replace and U on the right hand side by their renormalized values. The same reasoning 
applies with respect to this approximation as at the end of section 3, in particular for the 
important matrix elements at the Fermi surface. One finds 



(2) 

integral in (5.12) has a logarithmic divergence, which yields a logarithmic divergence of ^ if 



This formula is a very good approximation to (5.12) if e/j or are not too close to or e^ + U^. 
If both band energies become equal to or ef + U^, the approximate result diverges. In this 
special case the asymptotic behaviour of and U becomes important. One can show that the 

[2) 

■,q 

ejfc and approach or ef + U^. Such a divergence causes no problems since it is integrable 
as a function of the energy e^; also higher powers of V^"^^ remain integrable. 

(2) 

When calculating the Kondo temperature, we only need F^, ^ for the case where at least one 
of the band energies is equal (or at least very close) to the Fermi energy. Since or + 
can only be equal to the Fermi energy if the hybridization vanishes, we can use (5.13) in the 
following. At the Fermi surface this yields 

and with only one wave vector at the Fermi surface 

- ^..(0) n(0) U ^^^^^^^^ ^ ^^^^ ^ _ ^^^^^^^ ^^^^ . (5.15) 

Before proceeding with the calculation of the Kondo temperature in various cases, it is interesting 
to compare this result with the coupling obtained by the Schrieffer- Wolff unitary transformation 
in Ref. [9]. There one finds 

in particular at the Fermi surface 

<^. = V^M^j^£^y (5.17) 
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As a first remark we mention that for = both results are identical if one replaces the initial 
values of and in the result by Schrieffer and WolfT with the renormalized values and 
C/^. For gjt 7^ eq our result differs from the Schrieffer- Wolff result. 

The first problem with the induced spin-spin interaction in the formalism of Schrieffer and 
Wolff are the pole terms in vj^^^ ^ if or + lie in the conduction band. A second problem 
is apparent in the following limit 

This immediately implies that the effective band width in the corresponding Kondo problem is 
of order the conduction band width 

Deff OC D, (5.19) 

which is known to be wrong. It is quite obvious from Eq. (5.15) that both these problems do 
not show up in the flow equations approach. 

In order to obtain some more quantitative results, let us now discuss two particular regimes 
of the Anderson model. 



Symmetric Anderson model with 



< D 



In the symmetric case one has = and the relevant matrix elements of the spin-spin 

coupling are 

2 _ (UR)^ 



^ VkJO) Vg{0) ^. (5.20) 

[^1 ~) + [~^) 



This is depicted in Fig. 3 where it can be compared with the Schrieffer- Wolff result. For simplicity 
we have assumed Vq{0) = Vkp{0) for all wave vectors q in the diagram. Furthermore, we have 
replaced by in the Schrieffer- Wolff result, see also our discussion in section 6 for this 
point. The Kondo temperature depends mainly on the coupling at the Fermi surface (compare 
Eq. (5.6)) 

^-^fl-^V (5.21) 



kF,kF T^JJI \ J) J 

This is consistent with the Schrieffer- Wolff result except that we find an additional (usually 
small) correction term 2T/D. 

The main difference to the Schrieffer- Wolff result shows up in the effective band width of the 
associated Kondo problem. In the flow equations approach the effective band width is obviously 
proportional to the on-site interaction 

Deff OC (5.22) 



since the spin-spin coupling becomes ferromagnetic for \eq\ > U^/\/2 and decays to zero even 
further away from the Fermi surface. Since the spin-spin coupling induced by our unitary 
transformation is not constant (there is no physical reason why it should be), it is difficult to say 
quantitatively what the proportionality factor in Eq. (5.22) is: The Kondo problem is usually 
only treated for a constant spin-spin coupling with the conduction band electrons. Nevertheless, 
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we would like to estimate the proportionality constant in (5.22) for the special case of a constant 
density of states and constant ^^(0). This value can be compared with known results. To get 

(2) 

some rather approximate estimate one can e.g. replace ^ regarded as a function of by its 
value at the Fermi surface in an interval around the Fermi surface that has the same area as the 
original curve. This should give a lower bound on Dgg. We obtain 

DefF ~ 0.36?7^. (5.23) 
This result can be compared with the result from the Bethe-ansatz solution. There one obtains 

ksTK = -^eM-^iMeF)vll]k,))- (5-24) 
We identify the prefactor with the one in (5.6) and obtain 



Des = ^U = 0.443 U. (5.25) 

This is in good agreement with (5.23). Estimates similar to (5.23) can be obtained for an 
arbitrary density of states and arbitrary Vq{0). 



Valence fluctuation regime 

This regime shows characteristic new features of the asymmetric Anderson model. The spin- 
spin coupling is given by 

with the value at the Fermi surface 

MeF)viX = ^. (5.27) 
^^d 

Essentially the same remarks apply as in the previous section. The coupling is depicted in Fig. 4 
where it can be compared with the Schrieffer- Wolff result (again with Vg(0) = Vjt^(O) and with 
e^, replaced by ej, U^). As before the flow equations approach yields the correct scaling 
behaviour of the effective band width 

DeffOc|ef|. (5.28) 

It is important that the "renormalized" value of the impurity orbital energy enters in Eqs. (5.26) 
and (5.27). This behaviour is known from perturbative renormalization [10, 16] or numerical 
renormalization [6] which give the same value of the renormalized impurity orbital energy that 
we have calculated in Eq. (4.21). In this regime the Schrieffer-Wolff unitary transformation not 
only gives the wrong scaling behaviour of the effective band width D^e oc D, but also the wrong 
coupling at the Fermi surface since the initial value of the impurity orbital energy enters. 

At this point one can wonder about the temperature dependence of that has been found 
in Eq. (4.26). This effect seems to be unobserved in renormalization treatments. However, 
the maximum effect of non-zero temperature is to increase e^(T = 0) by a value of order F. 
According to Fig. 2 the smallest value of 
one might expect to see some influence of the temperature dependence of in this regime with 
F. But this is just the mixed valence regime (for ksT ^ T) where scaling breaks down 
anyway [16]. For this reason there is no contradiction. 



with < is = -F/tt, that is of order F. Thus 



.R 
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6 Comparison with other methods I (Schrieffer— Wolff unitary 
transformation) 

At this point we would like to explain in more detail why the Schrieffer-Wolff transformation 
and our continuous unitary transformation yield different results. At a first glance these two 
transformations are very similar. Our r] has the same general structure as the generator S in the 
Schrieffer-Wolff transformation H e^He^^ [9] with 



In both approaches the same more complicated interactions generated due to higher commutators 
are neglected. But as we have seen the two transformations still show important differences that 
have to be explained. 

First of all, let us mention that it is possible to construct a continuous unitary transformation 

(2) 

that reproduces the Schrieffer-Wolff result (5.16) for the coupling ^ with the initial values of 
erf and U replaced by their renormalized values. In fact one can achieve this by choosing the 
parametrization 

fieJ) = -l^ (6.2) 

instead of Eq. (4.8). Then the hybridization J{e,i) shows a very simple flow 

J(e,^) = J(e,0) exp(-/). (6.3) 

This is depicted in Fig. 5 where it can be compared with the behaviour for our choice of /(e,£) 
from Eq. (4.8). As compared to the original Schrieffer-Wolff transformation this is an improve- 
ment since now the renormalized parameters enter into the expressions for the induced spin-spin 

(2) 

interaction at the Fermi surface . In the valence fluctuation regime this is of importance 

as discussed in the previous section. 

Still the main differences between the Schrieffer-Wolff result and the flow equations are not 

(2) 

resolved so easily. That is the couplings ^ show the wrong high-energy cutoff and contain 
pole terms. Obviously our result for the spin-spin coupling diff'ers from the result by Schrieffer 
and Wolff already in second order in the hybridization V^. This is due to the fact that our 
transformation differs in this order from the Schrieffer-Wolff transformation. In principle it is 
of course possible to write our transformation in the form exp(5) too. S can be calculated from 
the expansion 

S= M r]{£) + - / di d£'[ri{£),r]{£')] + .... (6.4) 
Jo 2 Jo Jo 

In second order in the term containing the commutator of r) at two different values of the flow 
parameter becomes important. It is of the form 

k,q,a 

+ : C^q^^d'^_^dfjC]f^—fj : +C^^^C^ _^(i_o-C^(T ~ ^a^—a^q,—'^^k,a)- (6-5) 

Our choice of f{e,i) leads to a controlled expansion without any pole terms here: The worst 
divergencies occuring in the generated couplings are only logarithmic pole terms and are therefore 
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themselves integrable (compare the discussion in section 5). In contrast the parametrization (6.2) 
leads to non-integrable pole terms in (6.4) that cannot be interpreted as principal value integrals. 

Of course terms in order curing the divergencies could be introduced by hand into the 
original Schrieffer- Wolff generator S in Eq. (6.1). However, these would be divergent too, difficult 
to construct, and difficult to control in higher orders. In the flow equations formalism such terms 
appear naturally since r/ does not commute with itself for different values of i. Also the problems 
with pole terms are (as far as we have seen) resolved due to the introduction of a flow parameter 
that generates an asymptotic behaviour of the parameters and thereby "smears out" the pole 
terms. 



7 Comparison with other methods II (Anderson's "poor man's" 
scaling) 

Let us compare our method to perturbative renormalization in the spirit of Anderson's "poor 
man's" scaling approach [15]. We will briefly review the main features of this approach as applied 
to the Anderson impurity model by Haldane [16]. For simplicity we assume U ^ D ^ \eci\- The 
band width D is reduced to D — dD by perturbatively integrating out states with energies 
D — dD < \e\ < D. This results in the following scaling equations [16] 

dlnU Jo \ e + ea ^ - J 

'^^ 0{TJ\D)ID). (7.2) 



d InD 



In order to avoid confusion with the flow equations approach we have introduced the notation 

J^(e) =^ J(e, ^ = 0). The second equation leads to no nontrivial scaling behaviour. In the first 
equation the main contributions come from |e| D and Haldane finds 

'^'^ -- + 0{e,f{D)lD). (7.3) 



d \nD TT 

By scaling down io D k. ksT one obtains the effective impurity orbital energy Ej that is equal 



to our in (4.19). In the limit D ^ U ^ ej one argues that non-trivial scaling occurs only 
when D has been reduced to order U. In this case one obtains (4.22) for the effective impurity 
energy that enters into the induced spin-spin interaction. 

Now the scaling equation (7.1) can also be obtained in the fiow equations formalism. One 
has to choose 

(7.4) 

where is the initial band width. D{1) is a parameter (effective band width) that approaches 
zero monotonously as i — > oo. Instead of using the unfamiliar parameter i we can then try to 
work with the more familiar parameter D. Obviously (7.4) amounts to decoupling the high- 
energy modes by means of a unitary transformation as depicted in Fig. 5. It should come as 
no surprise that this yields flow equations in D reminiscent of the familiar scaling equations. 
Eq. (7.4) corresponds to a certain choice of f{e,i) but this will be of no importance for our 
discussion. 
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The flow equation (4.5) for in the Umit U ^ oo expressed as a function of D = D{£) takes 
the following form after a short calculation 

'^^'^ I' de{l+n{e))^—^j\e^) (7.5) 



d \nD J Ed — e de D 

to be integrated from D = to D = 0. For T = this obviously just reproduces Haldane's 
scaling equation (7.1). 

Now first of all this might make some approximations more plausible that we have used in the 
flow equations approach. In particular the restriction to a certain set of interactions is standard 
in a renormalization framework. It is well-known that the procedure of integrating out high- 
energy modes leads to additional interactions in the "poor man's" approach [16]. This will also 
include e.g. the induced spin-spin interaction. However, even these terms are not taken into 
account in the scaling approach.^ Instead one argues that on a low-energy scale the physical 
behaviour is determined by the scaling invariants like E*^. 

Though intuitively quite appealing, the parametrization (7.4) leads to immediate problems 
in the flow equations formalism since we do not use similar approximations. The first problem 
is already apparent from (7.2): Why should the infinitesimal unitary transformations stop when 
D is of order ksTl Such a cutoff is necessary, otherwise Eq. (7.2) leads to divergent behaviour 
as D ^ 0. 

An even worse problem shows up when comparing the self-consistency equations (4.24) and 
(4.25) with the actual solution of the differential equations in D. Let us for simplicity consider 
the symmetric Anderson model with (compare Eq. (4.13)) 



However, the solution of the differential equation in D leads to U{1 = oo) = as can be shown 
rigorously for the semi-circle case. Intuitively what happens is that the modes with energies 
larger than ^7/2 that are being integrated out push U closer to the Fermi surface. Before the 
low-lying modes can push U away again, one has reached [/ = 0. Just this behaviour is avoided 
by the parametrization (4.7), compare also Fig. 5. This is the reason why the decoupling in our 
case can be performed all the way down to £ = oo without introducing a condition like D « ks T. 
The self-consistency conditions (4.24) and (4.25) have been obtained from (4.10) and (4.11) by 
taking ^ = 0. Thus the self-consistency conditions fail if the principal value integrals in (4.10) 
and (4.11) are not well-defined for some value of I. This happens with the parametrization 
(7.4). For I such that D{1) = or D{i) = + the integrals diverge, a familiar band-edge 
problem. 

The reason why this phenomenon is unobserved in the scaling approach is that it is due 
to neglected terms like on the right hand side of Eq. (7.2). But if one attempts to study the 
behaviour when U becomes of order D, these terms cannot be neglected. Another way of saying 
this is that the scaling equations do not reproduce a standard perturbational expansion like e.g. 
Eq. (7.6) and a behaviour like Eq. (4.19) within one framework. The flow equations manage this, 
but not with a parametrization like (7.4). One might hope that they can provide a tool to study 
crossover phenomena. 



^In this sense our approximations are clearly better though less familiar. 
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8 Conclusions 



We have shown that using a continuous unitary transformation it is possible to eliminate the 
hybridization term in the single impurity Anderson model. Similar to the well-known Schrieffer- 
Wolff transformation, an additional spin-spin interaction and pseudo-spin interaction are gener- 
ated. Furthermore the impurity energy and the on-site interaction are renormalized. The main 
approximation was to neglect additional interactions generated by the transformation if the spin- 
spin interaction is included in the Hamiltonian. Instead we simply integrated the coefficient of 
this interaction to obtain the new coupling constants. Some justification for this procedure is 
that by introducing the spin-spin interaction in the Hamiltonian right from the beginning, we 
do not obtain additional contributions to the impurity energy and the interaction U. Such con- 
tributions occur only if even more complicated additional interactions are taken into account. In 
this manner we have restricted ourselves to some minimal but "consistent" set of interactions. 

This procedure yields coupled non-linear flow equations for the impurity energy and the in- 
teraction. It turns out that these equations can be approximately solved by self-consistency 
conditions for the renormalized values of the impurity energy and the interaction. This ap- 
proximation can be controlled by comparing the result with a numerical integration of the flow 
equations and it turns out that the approximation is very good. It becomes exact in the limit 
of vanishing interaction U. We discussed several regimes in the parameter space for a special 
coupling J(e, 0) = -^^/D"^ — e^. We do not expect that the general features depend on the 
details of J(e, 0), but only on F and D. The method works in the whole parameter space as long 
as the hybridization is not the largest energy scale in the system. Although we mainly discussed 
the case where J(e, 0) is a semi-circle, the method can be applied to any other J(e, 0) as well. 

The Kondo regime of the Anderson impurity model has been investigated by mapping it 
onto the Kondo Hamiltonian with a unitary transformation. In contrast to the Schrieffer-Wolff 
unitary transformation our approach of infinitesimal unitary transformations generates all the 
correct parameters of this Kondo Hamiltonian. In particular no pole terms appear in the spin- 
spin coupling. An important point is also that the correct high-energy cutoff enters in all the 
expressions for the renormalized or induced interactions. One should expect that such a unitary 
transformation exists since after all the Anderson impurity model has been introduced to study 
dilute magnetic alloys which exhibit the Kondo effect. 

Quite generally the flow equations seem to provide a good framework to investigate such 
induced interactions like the spin-spin interaction here. Another interesting problem of this 
type is the problem of electrons interacting with phonons in a solid. In this case the well- 
known Frohlich transformation decouples the electrons from the phonons [17]. It generates an 
additional effective interaction between the electrons. In some region of the parameter space 
this interaction is attractive and is responsible for the formation of Cooper pairs. But again 
the Frohlich transformation is ill-defined near resonances. Divergencies as in the case of the 
Schrieffer-Wolff transformation occur. This problem has recently been studied by Lenz [18] 
using continuous unitary transformations. As in our case his transformation is well-defined and 
the attractive interaction between the electrons does not show the usual divergencies. 

The examples discussed so far show that the method of continuous unitary transformation has 
a wide range of applications. They provide a useful tool to simplify a given Hamiltonian. In the 
examples discussed so far [1, 7, 8, 18] a given Hamiltonian was transformed into a block-diagonal 
form by eliminating some of the coupling constants. This generates additional couplings and the 
original problem is mapped thereby to a different problem which can be investigated further. 
The main problem is that one has to neglect some of the additional couplings to obtain a closed 
set of equations. Up to now we are not able to present a general method that allows to estimate 
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the effect of this approximation. The reasoning should probably be similar to the discussion of 
irrelevant terms in renormalization theory. But the comparison of our results with well known 
results obtained by other methods show that our approximation works well. In particular we can 
treat a strongly correlated system like the Anderson impurity model, where the natural small 
parameter is the hybridization F and the rest of the Hamiltonian is not of simple quadratic form 
but contains the large interaction U. 

A partially open problem is to calculate expectation values. To do this one has to apply the 
continuous unitary transformation to the observable as well. If one is only interested in static 
properties, this seems to be no problem. One has to solve the flow equations for the observable, 
which needs additional approximations for the observable similar to the approximations made 
for the Hamiltonian itself. Since we neglect only normal ordered terms, which do not contribute 
to the expectation value, this should be possible [1]. But if one wants to calculate dynami- 
cal properties like time-dependent correlation functions, it is yet not clear whether the same 
approximations can be applied. This problem has been discussed to some extent in Ref. [8]. 

One situation where one needs to calculate dynamical properties is the Hubbard model in 
the limit of infinite space dimensions [19]. It can be mapped onto a single impurity Anderson 
model with an additional self-consistency condition for the local density of states [20]. To apply 
our method to this problem, we have to calculate the single-particle Green's function, which is 
a dynamical problem. Also it would be important to attempt to study the crossover behaviour 
in the Kondo problem for low temperatures using our approach. Work in these directions is in 
progress. Other interesting problems are for example the single impurity Anderson models with 
two or more channels. 
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Figure captions 

Fig. 1. Various regimes for the solution of the U = Anderson impurity model. 

Fig. 2. The behaviour of the renormalized impurity orbital energy close to the Fermi surface as ob- 
tained from Eq. (4.19). The dashed line corresponds to the solution of the self-consistency 
condition, the full line is deduced from the numerical solution of the differential equations. 

Fig. 3. Induced spin-spin interaction in the symmetric case. 

Fig. 4. Induced spin-spin interaction in the [/ — )■ oo limit. 

Fig. 5. Sketch of J(e, I) for the symmetric Anderson impurity model for different values of ^, 

£2> ii> io = or D2 < Di < Do = . J(e, 0) is the semi-circle. 

In a) J{e,£) is shown for the improved Schrieffer- Wolff transformation introduced in (6.3). 

b) corresponds to the poor man's scaling approach in the context of the flow equations, 
compare sect. 7. 

In c) the parametrization of J(e, I) as chosen in this paper is depicted. The flow of U {£) 
has been neglected in the diagram, in fact this would lead to an algebraic decay of J(e, i) 
also for e = ±U^/2. 
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